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Abstract. We study the surfaee photometry for 13 
single-nucleus ultraluminous IRAS galaxies (ULIRGs), se- 
lected from a parent sample of 58 galaxies obtained in 
a Hubble Space Telescope snapshot survey. We find that 
these galaxies can be classified into three classes according 
to their surface photometry. The surface brightness pro- 
files of the four objects in the first class are well fitted by 
the R^ law. The isophotes are all disky at i? < lh~^ kpc , 
consistent with the molecular disks/rings found in nearby 
ULIRGs from CO observations. Each of the four galax- 
ies in the second class has an inner R^ component plus 
an outer extension. Remarkably all these four galaxies 
are Seyfert 1 galaxies with luminosities in the quasar 
regime and with relatively narrow permitted and strong 
Fell emission lines. The remaining five galaxies fall into 
the third class; these objects have surface brightness pro- 
files that deviate to various degrees from the R* law; in- 
deed, one is fitted perfectly by an exponential law. We also 
present new spectroscopic observations for some of these 
galaxies. We establish the redshift of IR 09427-1-1929 to 
be 0.284, instead of 0.149 as adopted in the literature. 
These observations firmly support the idea that merging 
of disk galaxies produce ellipticals. Our data also strongly 
suggest that the formation of QSOs may be an integral 
part of elliptical galaxy formation (at low redshifts). 

Key words: Galaxies: photometry Galaxies: interactions 
Galaxies: elliptical Galaxies: Seyfert Galaxies: nuclei 



Send offprint requests to: Z. Zheng (zhengz@bac.pku.edu.cn) 
^ BAG is jointly sponsored by the Chinese Academy of Sciences 
and Peking University. 

* This research was based on observations obtained with the 
NASA/ESA Hubble Space Telescope through the Space Tele- 
scope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 



1. Introduction 

Ultraluminous IRAS galaxies (ULIRGs) , characterized by 
their high far-infrared luminosities Ljr > IO^^Lq (for 
Hq — 50 km s~^ Mpc~^), are the most luminous galax- 
ies in the local universe (see Sanders & Mirabel 1996 
for a review). Most ULIRGs show signs of interaction 
and/or merging. These processes are believed to trigger 
starbursts and central nuclear activities which in turn pro- 
duce intense infrared emissions via dust absorption and 
re-emission. Some ULIRGs do not show prominent inter- 
action signatures (at least from ground-based imaging), 
but these could be advanced mergers, as is the case for the 
ULIRGs with highest luminosities (Lir, > 2.25 x IO^^Lq, 
for Ho = 50 km s'^ Mpc^^). These ULIRGs may repre- 
sent an important stage in formation of QSOs and power- 
ful radio galaxies, and may also present a primary stage in 
the formation of elliptical galaxies. (Sanders et al 1988a; 
Sanders & Mirabel 1996; Melnick & Mirabel 1990). 

Although the ultimate fate of the ULIRGs is not yet 
fully understood, more and more observations show that 
the remnants of the mergers resemble elliptical galaxies 
or SO galaxies. For example, the surface brightness profile 
of IRAS 20551-4250 is reasonably well fitted by the Ri 
law (Johansson 1991). Arp 220 and NGC 6240, two near- 
est ULIRGs, both have densities, velocity dispersions and 
central surface brightness distributions consistent with the 
fundamental plane of elliptical galaxies (Doyon et al 1994). 
Baker & Clements (1997) investigated the old stellar pop- 
ulation for 10 nearby ULIRGs based on deep near-infrared 
imaging and concluded that 8 out of 10 ULIRGs show 
signs of elliptical structure. 

On the theoretical side, Toomre & Toomre (1972) first 
successfully reproduced the observed bridges and tails us- 
ing restricted three-body simulations of two colliding disk 
galaxies. In 1977, Toomre put forward the hypothesis that 
all elliptical galaxies could be remnants of merged disk 
galaxies. Later simulations improved the spatial resolu- 
tions, explored more phase space of the encounter geome- 
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try and incorporated more realistic treatments of the pro- 
genitors, adding components such as bulges, massive dark 
matter haloes and gas; crude treatment of star formation is 
also attempted (e.g. White 1978; Farouki & Shapiro 1982; 
Negroponte & White 1983; Barnes 1992; Hernquist 1992; 
Heyl et al 1994; Barnes & Hernquist 1996; Weil & Hern- 
quist 1996; Walker et al 1996; Mihos & Hernquist 1994). 
These simulations conclude that a wide range of initial 
conditions can result in the formation of early-type-like 
galaxies (see de Zeeuw & Franx 1991 for a review). Besides 
the global elliptical-like characteristics of the remnants of 
simulated mergers (e.g. the i?* law of the surface bright- 
ness profile), some simulations also attempt to reveal the 
fine structure of the remnants such as the shape of the 
isophote, loops, shells and tails (see e.g. Heyl et al 1994; 
Springel & White 1998). As mergers of galaxy pairs may 
not account for all observed properties of ellipticals, mul- 
tiple mergers have been considered (Mamon 1987, Barnes 
1984, Barnes 1985, Barnes 1989, Schweizer 1989; Weil & 
Hernquist 1996); multiple merging may occur in compact 
groups such as those found by Hickson (1982, 1993). 

While the simulations indicate that merging of two 
or more disk galaxies can easily lead to the formation 
of ellipticals, it is clearly important to test these pre- 
dictions empirically. For this purpose, ULIRGs are ideal 
since they are closely related to galaxy merging and in- 
teractions (e.g.. Murphy et al 1996, Clements et al 1996). 
In this paper, we analyze the surface photometry for 13 
ULIRGs observed with the Hubble Space Telescope (HST) 
Wide Field Planetary Camera 2 (WFPC2), and examine 
whether their profiles resemble the R^ law, as one would 
expect if they ultimately form elliptical galaxies. The res- 
olution of HST provides us a powerful way to study the in- 
ner profiles and the fine structure details in these galaxies. 
We will also address the important question whether the 
surface photometry of these galaxies are associated with 
the central active galactic nuclei (AGN) phenomenon. The 
structure of the paper is as follows. In section 2, we de- 
scribe how we select our ULIRG sample. Data reduction 
including surface photometry is described in Section 3. In 
Section 4, we classify these galaxies according to their sur- 
face brightness profiles. In Section 5, we discuss the spec- 
troscopic properties of these galaxies, using both newly 
obtained spectra and data from the literature. In Section 
6, we summarize and discuss the implications of our re- 
sults. Throughout this paper we adopt an Einstein-de Sit- 
ter universe (Jig = 1) Sind denote the Hubble constant as 
Hq = lOO/i km s^^ Mpc^^. Since our objects are at rela- 
tively low redshift (z ^ 0.35), the assumption about the 
density parameter (rig) is not critical, and 1" corresponds 
to roughly Ih^^ kpc (z/0.1) for z <, 0.4. 

2. Sample Selection 

In this study, we use archive images from an HST snapshot 
survey of ultraluminous IRAS galaxies (Borne et al 1998). 



This snapshot survey was taken using the WFPC2 in the 
I bandpass (F814W). The objects observed were mainly 
selected from the following samples: the bright samples of 
Sanders et al (1988a,b) and Melnick & Mirabel (1990), and 
the QDOT sample (Leech et al 1994; Lawrence et al 1999). 
For each target in this survey, two 400s exposures were 
taken. Most of target galaxies are centered on the Wide 
Field Camera chip 3 (800x800 pixels, 0.0996" per pixel) 
and a few in the Planetary Camera chip (800x800 pixels, 
0.0455" per pixel). For details of this survey, see Borne et 
al (1998). In total, there are 58 images available to us from 
March 1996 to Jan. 1997 (except that IR 10026-^4347 was 
taken in May 1997). The redshifts of these objects range 
from 0.04 to 0.35. 

Table 1. Journal of Observations 



Target Name 


RA(2000) 


Dec(2000) 


Date 


IR 09427-^1929 


09:45:29.1 


19:15:50.0 


23/05/96 


IR 17432-5157 


17:47:09.9 


-51:58:44.0 


27/10/96 


IR 23140-^0348 


23:16:35.2 


04:05:17.0 


03/09/96 


IR 04025-8303 


03:57:11.3 


-82:55:16.0 


02/07/96 


IR 00509-H1225 


00:53:34.8 


12:41:36.1 


29/08/96 


IR 00276-2859 


00:30:04.1 


-28:42:25.6 


05/10/96 


IR 020554-0835 


02:08:06.7 


08:50:03.6 


05/10/96 


IR 100264-4347 


10:05:41.9 


43:32:39.4 


01/05/97 


IR 14378-3651 


14:40:59.4 


-37:04:33.2 


02/09/96 


IR 105594-3845 


10:58:39.3 


38:29:06.4 


19/06/96 


IR 20176-4756 


20:21:11.1 


-47:47:07.1 


28/09/96 


IR 09039-H0503 


09:06:34.1 


04:51:28.3 


17/11/96 


IR 23242-0357 


23:26:50.2 


-03:41:05.5 


24/08/96 



The near infrared observation (F814W) minimizes the 
effect of the absorption by gas and dust in the target galax- 
ies because the dust extinction at longer wavelengths is 
smaller. Moreover, near infrared imaging is not sensitive 
to young stellar populations. The I-band imaging is there- 
fore more suitable for our primary goal of understanding 
how the old stellar populations relax dynamically during 
the merging process. 

The ULIRGs in the Borne et al (1998) sample exhibit 
a variety of morphologies. Some have only a single nucleus 
while others have two or multiple nuclei and display com- 
plex and disturbed configurations. These morphologies are 
likely connected to different stages of merging and inter- 
action. In this paper, we will focus on studying properties 
of the late-products of merging and relaxation, therefore 
we choose a subsample of galaxies that are dominated by 
a single nucleus; these single-nucleus galaxies are selected 
by checking the contour plot for each source. The tar- 
gets in our sample all show merging features, such as tidal 
tails, plumes, shells and other fine structures. The galax- 
ies selected through this way are probably all mergers at 
the late stage of interacting/merging (§4, see also Wu et 



al 1998). In total, 13 targets are chosen out of the 58 im- 
ages. The journal of observations for these targets is listed 
in Table 1. 



3. Data Reduction 

All observations were preprocessed "on the fly" through 
the standard Space Telescope Science Institute pipeline as 
described by Holtzman et al (1995). The images were cal- 
ibrated (including bias subtraction and flat fielding) with 
the most up-to-date version of the routine reference files 
provided by the Institute at the time the images were 
taken. After correction of hot pixels, our image processing 
mainly include cosmic-ray removal and sky subtraction, 
which we describe below. 



3.1. Cosmic-ray Cleaning 

With its high altitude, HST is free of atmospheric seeing 
effects, however this also means that it is much more sus- 
ceptible to cosmic ray hits than ground-based telescopes. 
So it is essential to carefully remove the cosmic rays. 

Since for each of our object there are two 400s ex- 
posures that are in general well aligned, the cosmic-ray 
removal is relatively straightforward. We use the stan- 
dard IRAF task "combine" with a (conventional) 3a clip- 
ping threshold to detect the cosmic-rays. The cosmic-ray 
free images are then co-added. The combined image has a 
higher signal-to-noise ratio than each single image, usually 
by a factor of v^ except for some individual pixels. 

3.2. Sky Subtraction 

The background of each combined image was fitted by per- 
forming the IRAF task "imsurfit" . Since the field of view 
for the Wide Field Camera is small (about 80" x 80"), 
it is reasonable that the sky background has a relatively 
small variation across the field. After several trials, we 
found that it is sufRcicnt to fit the sky background using 
a two-dimensional second-order Legendrc function, i.e., a 
"declined plane" . Some rectangular regions (usually each 
includes more than 2000 pixels) were chosen from the com- 
bined image to sample the sky background. The median 
value of each region was used for the fitting. 

The sky-subtracted image is then ready for photomet- 
ric purposes. 

3.3. Surface Photometry 

We perform the surface photometry using the ISOPHOTE 
package in STSDAS in the IRAF environment. 

Foreground stars and diffraction spikes caused by 
bright stars or central nuclei in each image were care- 
fully masked out to eliminate their effect on the photom- 
etry. We use the SExtractor galaxy photometry package 
(vl.2bl0b) (Bertin & Arnouts 1996) to help identify the 



stars. A circular region and a triangular region are used to 
mask a star and a spike, respectively. Pixels inside these 
regions are given zero weighting in the ellipse fitting pro- 
cedure which we describe below. 



3.3.1. Isophotc Ellipse Fitting 

For all galaxies, the IRAF task "ellipse" was used to derive 
the isophotal morphological parameters. 

The "ellipse" task fits ellipses which best reproduce 
the observed isophotes of an image (Jedrzejewski 1987). 
The best parameters are found by an iterative procedure. 
In each run of the isophotc fitting, a number of parame- 
ters are used to describe the fitted ellipse. These include 
the ellipse center (xQ^yo), semi-major axis a, ellipticity 
e (e = 1 — b/a, where b is the semi- minor axis), posi- 
tion angle cf) (the azimuthal angle of the major axis, mea- 
sured counter-clockwise from the +y direction in the im- 
age here) . To describe the deviation of isophotes from per- 
fect ellipses, the difference between the isophotal radius 
and the best ellipse fit is usually expanded in Fourier se- 
ries. The coefficient of cos 40 term in the expansion, B4, 
is given particular importance in describing the deviation 
(e.g. Carter 1977, Jedrzejewski 1987; Kormendy & Ben- 
der 1996): if B4 is positive, the isophotc is more elongated 
along the major axis than the best fitting ellipse and is 
said to be disky; on the other hand, if B4 is negative, 
the isophote will appear rectangular (boxy). The inten- 
sity around each isophote is measured by averaging the 
fiux azimuthally. The ellipse center, ellipticity and posi- 
tion angle are allowed to change in this work for two rea- 
sons: First, ellipticity changes and isophote twisting (i.e. 
changes of position angle) have already been noticed in 
elliptical galaxies. Second, our sample mainly includes in- 
teracting/merging galaxies. Behavior of the changes in the 
parameters may reveal interesting structure details. 

3.3.2. Fitting the Surface Brightness Profile 

In order to address the questions mentioned in the intro- 
duction, we checked the one-dimensional surface bright- 
ness profiles obtained through isophote ellipse fitting. 

We fit the radial profiles with a Scrsic (1968) law (Mar- 
Icau & Simard 1998): 

I{r) ^ h exp(-fe[(^)i/" - 1]), b ^ 1.9992n - 0.3271, (1) 

Re 

where I{R) is the surface brightness at radius R, and Re 
is the effective radius within which half of the light is en- 
closed (Capaccioli 1989). The de Vaucouleurs R^ profile, 
which is usually used to describe elliptical galaxies and 
bulges, is a special case of (1) with n = 4, and the expo- 
nential profile, which is used to describe galactic disks, is 
also a special case with n = 1. 
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Table 2. Properties for the Ultraluminous IRAS Galaxies 



Target Name 


Redshift 


Mi-5 log h 


Log(ft^-L/ir/L0) 


Re(ft"^kpc) 


Spectral Type 


Class 


Comments 


IR 09427+1929 


0.28400 


<-23.21 


12.01 


2.75 


Syl° 


1 


Strong Fell 


IR 17432-5157 


0.17500 


-21.11 


11.73 


4.99 


LINER" 


1 




IR 23140-F0348 


0.21980 


-22.94 


11.79 


2.47 


LINER° 


1 


Radio Galaxy 


IR 04025-8303 


0.13660 


<-22.41 


11.55 


0.94 


Syl' 


1 




IR 00509-^1225 


0.06114 


<-22.72 


11.21 


2.69 


Syl" 


2 


Strong Fell 


IR 00276-2859 


0.28000 


<-23.70 


12.04 


3.25 


Syl" 


2 


Strong Fell 


IR 02055+0835 


0.34500 


<-23.31 


12.37 


0.80 


Syl' 


2 


Strong Fell 


IR 10026+4347 


0.17800 


<-23.11 


11.59 


0.84 


Syl' 


2 


Strong Fell 


IR 14378-3651 


0.06760 


-21.17 


11.72 


- 


LINER*^ 


3 




IR 10559+3845 


0.20660 


-22.32 


11.74 


- 


HIF 


3 




IR 20176-4756 


0.17810 


-21.51 


11.81 


- 


not AGN' 


3 




IR 09039+0503 


0.12500 


-21.80 


11.70 


- 


LINER° 


3 




IR 23242-0357 


0.18900 


-21.64 


11.52 


- 


HIF 


3 





NED; '' Spectra of QDOT; ' Lawrence et al 1999; *= Kim et al 1998; ° Our observation 



4. Results 

Table 2 gives the properties for these galaxies, includ- 
ing their redshifts, magnitudes, infrared luminosities and 
spectral classifications. The integrated apparent WFPC2 
magnitude is obtained from aperture photometry in the 
sky-subtracted and star-masked image using the formula 
from Holtzman et al (1995); the apparent magnitude is 
then converted into absolute magnitude using luminosity 
distances for our cosmology (no k-correction has been ap- 
plied). In the case of saturated images, only upper limits 
(in absolute magnitude) are given. Figs. 1-13 (left panels) 
show the HST images for all thirteen galaxies. Each im- 
age has 750 X 750 pixels, corresponding to 74.5" on a side. 
The middle panels show the contour levels for a smaller 
region centered on each target, highlighting its environ- 
ments. The right panels show the variations of the surface 
brightness, ellipticity, position angle and B/^/a as a func- 
tion of i?4 . From the radial surface brightness profiles we 
classify these galaxies into three classes: 



I. The first class includes 4 galaxies; the surface bright- 
ness of each galaxy in this class is well fitted by the R^ 
law with some small deviations at ^ h^^kpc. These 
four galaxies have round appearances. 
II. The second class includes 4 galaxies. The surface 
brightness profile for each galaxy in this class has two 
components: an inner R^ component and an outer 
component. There is a bright central nucleus in each 
of these galaxies. The best-fit effective radii (for the 
i?4 component) for class I and II galaxies are given in 
Table 2. 
III. The remaining 5 galaxies fall into the third class. These 
galaxies have surface brightness profiles that deviate 
significantly from the R^ law. 



In the following, we will discuss the properties of these 
three classes in detail. 



4-1. Class I 

Four galaxies (IR 09427+1929, IR 17432-5157, IR 
23140+0348, IR 04025-8303) belong to class I. The 
surface brightness profiles for these galaxies are well de- 
scribed by the R^ law, from the very inner part to 
the outer part as shown in Fig. 1 to Fig. 4. For IR 
17432—5157, the fit is reasonable from 0.5 /i~^kpcto 
9.2h-^ kpc ; for IR 23140+0348, the fit is excellent out to 
13/i^^ kpc . Another common characteristic of these galax- 
ies is that they have round appearance. For each galaxy, 
the maximum ellipticity of the isophotes is less than 0.30. 
The cllipticities at the effective radius are even smaller. 
They are only 0.2, 0.1, 0.1, and 0.1 for IR 09427+1929, 
IR 17432-5157, IR 23140+0348, IR 04025-8303, respec- 
tively. 

However, there are slight deviations from the exact 
Ri law. Small excesses over the best-fitted curve can be 
seen in all surface brightness profiles (right panels in Figs. 
1-4) at ~ (1 — 3) ft." ^ kpc. The bump position in each 
galaxy coincides with the maximum of the ellipticity e. 
This coincidence reflects small structures in these galax- 
ies. Furthermore, it is very interesting that all four galaxies 
have disky isophotes {B4 > 0) in the very inner regions 
(< l/i~^ kpc), and the position angle change is very large, 
(^ 100° for some galaxies). CaroUo et al (1997) analysed 
HST WFPC2 images for 15 elliptical galaxies. In both 
F555W and F814W filters, the position angle changes are 
usually less than 20°. Therefore it appears that although 
the light profile is described by the i?* law, the ellipses 
still show substantial twists. 



IR 09427+1929 



EPS files of the images are available at 
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Fig. 1. The left panel shows the HST Wide Field Camera image for IR 09427+1929 (the object labeled as GO). The 
image shown is 74.5" on a side. Three additional galaxies (Gl, G2 and G3) are within the field. The ellipse shows 
the error box of the IRAS position (cf. §5.1). The middle panel shows the contours centered on IR 09427+1929; the 
contour levels are 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024 ADU, respectively. The thick horizontal bar indicates a scale 
of 5h~^ kpc . The north and east directions are indicated at the top right in the middle panel. The right panels show 
the variations of the ellipticity, position angle, surface brightness, and ^4/0 (see §3.3) as a function of R^. The best 
Ri fit to the surface brightness profile is given by the straight line. 



IR 17432-5157 



EPS files of tfie images are avaiiabie at 
http://vega.bac.pltu.edu.cn/~zhengz/preprint.html 
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Fig. 2. Same as in Fig. 1, for IR 17432—5157. At least one additional galaxy are within the field. Due to its low 
Galactic latitude (—12°), many stars are also present. The contour levels are 1.9, 3.8, 7.6, 15.2, 30.4, 60.8, 121.6, 243.2, 
486.4 ADU, respectively. 



From Figs. 1-4, it can be seen that all four galaxies ex- 
hibit signatures of merging such as short tails, plumes in 
their outer parts. Except for IR 23140+0348 , all the other 
galaxies have some other galaxies within the same field of 
view, but these could be due to chance alignment (see Sec- 
tion 5). The round appearance of these galaxies provides 
a clue to the merging environment. Numerical simulations 
(e.g. Weil & Hernquist 1996) found that remnants of mul- 
tiple mergers are more round than the mergers from pair 
stellar disks (when projected onto the sky). So the round 
appearance of the class I galaxies suggests that they may 
be remnants of multiple mergers. Such a possibility is sup- 
ported by the high percentage of multiple nuclei galaxies 



in our parent sample of 58 galaxies. About 50% of these 
have more than three nuclei. Furthermore, for some of the 
double- nucleus galaxies, each of the two merging nuclei 
could contain a pair of galactic nuclei in high-resolution 
images, similar to Arp 220 (Taniguchi & Shioya 1998). 
Therefore, the percentage of the multiple-nucleus galaxies 
is likely higher than 50% for our sample, supporting the 
notion that ULIRGs may preferentially reside in groups of 
galaxies (Wu et al 1998). 

The physical significance of the disky isophotes at 
the inner radius (< 1/i^^kpc) and the excesses at (~ 
3/i~^kpc) is not clear. Multiple-waveband observations 
of ULIRGs reveal the presence of massive starbursts in 
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IR 23140+0348 



EPS files of the images are available at 
http;//vega. bac.pku.edu.cn/~zhengz/preprint. html 
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Fig. 3. Same as in Fig. 1, for IR 23140+0348 (labeled as GO). Faint features are visible around GO. The contour levels 
are 1.64, 3.28, 6.56, 13.12, 26.24, 52.48, 104.96, 209.92, 419.84, 839.68, 1679.36 ADU, respectively. 



IR 04025-8303 



EPS files of the images are available at 
http://vega.bac.pku.edu.cn/~zliengz/preprint.html 





(R/h"'kpc) 



Fig. 4. Same as in Fig. 1, for IR 04025—8303. At least one additional galaxy (Gl) is within the field. The contour 
levels are 1.55, 3.1, 6.2, 12.4, 24.8, 35, 49.6, 70, 99.2, 150, 198.4, 396.8, 793.6, 1587.2, 3174.4 ADU, respectively. 



a molecular ring or disk at several hundred pc to kpc in 
the circumnuclear regions (e.g., Downes & Solomon 1998). 
The positions of the inner disky isophotes may reflect the 
presence of such small disks at the center. Such nuclear 
disks are also found in numerical simulations where gas is 
driven inward to the central kpc region and cool efficiently 
at early stage of galactic encounters (Icke 1985; Noguchi 
1988). Therefore, the inner "disky" structure could be the 
relics of nucleus starburst rings or disks. The significance 
of bumps at ~ 3h~^ kpc is less clear, since even "normal" 
elliptical galaxies do not follow the R^ law perfectly (e.g., 
Saglia et al 1997). These departures could be produced 
in the merging process by tidal force and/or dynamical 
friction. 



4.2. Class II 

Figs. 5-8 show the four galaxies in class II. Surface bright- 
ness profiles of all these galaxies can be fitted by the R^ 



law from ~ Ih^^ kpc out to ~ 5h^^ kpc. The outer 
parts of these galaxies show clear excess over the R^ law. 
We have tried a two-component fitting with an inner i?* 
component and an outer exponential component. We find, 
however, that the fitting is unsatisfactory, meaning that 
the outer extension is more complex than a simple expo- 
nential. Disky isophotes are found in the inner region of 
each galaxy. However, this is clearly an artefact caused 
by the saturated central nuclei, since the diagonals of the 
quasi-diamond shaped isophotes are almost along the di- 
agonals of the images. 

Although our classification of these galaxies are based 
purely on their surface brightness profiles, remarkably, we 
found that all the four galaxies host luminous bright nu- 
clei and are saturated in the frames. All the galaxies were 
classified as Seyfert 1 galaxies spectroscopically (Lawrence 
et al. 1999; Vader & Simon 1987; Moran et al 1996), and 
all have luminosities in the QSO regime. For details, see 
Sect. 5.2. 



IR 00509+1225 



EPS files of the images are available at 



http://vega.bac.pku.edu.cn/~zhengz/preprint.html 




.'„nH»-"'ifini;l 'mi 




N 

- 1/n=0.25 ^ 





1 1.2 1.4 I.e I.B 1 l.g 1.4 1.6 1.B 

Arcseconds , 

(R/h"'kpc)"* 

Fig. 5. The left panel shows the HST Wide Field Camera image for IR 00509+1225 (the object labeled as GO). The 
image shown is 74.5" on a side. A satellite galaxy (Gl) and a foreground star SI are visible. The middle panel shows 
the contours centered on IR 00509+1225; the contour levels are 1.7, 3.4, 6.8, 13.6, 27.2, 54.4, 108.8, 217.6, 435.2, 870.4, 
1740.8 ADU, respectively. The thick horizontal bar indicates a scale of 5h~^ kpc . The north and east directions are 
indicated at the top right in the middle panel. The right panels show the variations of the ellipticity, position angle, 
surface brightness, and B^/a (see §3.3) as a function of R^ . The best R* fit to the surface brightness profile is given 
by the straight line. 

IR 00276-2859 



EPS files of the images are available at 
http://vega.bac.pku.edu,cn/~zhengz/preprint.html 
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Fig. 6. Same as in Fig. 5, for IR 00276—2859. A number of galaxies (G1-G5) and a star (SI) are within the field. The 
contour levels are 1.65, 3.3, 6.6, 13.2, 26.4, 52.8, 105.6, 311.2, 622.4, 1244.8, 2489.6 ADU, respectively. 



IR 00509+1225 It has a small galaxy with a projected 
distance of 13/i^^kpc away to the west. Two armlike 
features are clearly visible. The optical feature from 
ground based deep observation looks like that of Mrk 
231 (for details see Surace et al 1998). These armlike 
features could be tidal tails resulting from merging. 

IR 00276—2859 There are a number of small galaxies 
within the same field. In the contour plot, a small 
galaxy (Gl) to the southeast can be detected. Gl 
causes the outer contours of IR 00276—2859 to pro- 
trude in its direction, indicating that these two galaxies 
are physically interacting. 

IR 02055+0835 This galaxy has an extension to the 
west. Two objects Gl and G2 are within the same 



field. Both objects are within a projected distance of 
20h~^ kpc and could be physically associated with the 
IRAS galaxy. Many small and fainter objects are also 
visible within 20". 
IR 10026+4347 This galaxy is an IR QSO and is po- 
tentially a member of one important class of transi- 
tion objects. Detailed discussions for this object will 
be found in Xia et al (1999). 



4.3. Class III 

The five galaxies in this class have different surface bright- 
ness profiles, and so we discuss them in turn. 
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IR 03055+0835 



EPS files of the images are available at 5 _ 

http://vega.bac.pku.edu.cn/~zhengz/preprint.html 
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Fig. 7. Same as in Fig. 5, for IR 02055+0835. The contour levels are 1.7, 3.4, 6.8, 13.6, 27.2, 54.4, 108.8, 217.6, 435.2, 
870.4 ADU, respectively. Two nearby objects (Gl and G2) are clearly visible. 



IR 10026+4347 



EPS files of tfie images are available at 
http://vega.bac.pku.edu.cn/~zhengz/preprint.html 
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Fig. 8. Same as in Fig. 5, for IR 10026+4347. The contour levels are 1.75, 3.5, 7, 14, 28, 56, 112, 224, 448, 896, 1792 
ADU, respectively. 



IR 14378—3651 It is a galaxy with an asymmetric ex- 
tension, being more diffuse to the northwest, as can be 
seen in Fig. 9. A small galaxy (Gl) is visible about 
4/i~^ kpc (projected) to the southwest. This galaxy 
is still within the envelope of IR 14378— 3651 and 
could be physically associated with IR 14378—3651 . A 
ripple-like material can be separated from the south- 
east part of the extension. A loop structure can be 
seen in the inner part which is similar to spiral arms. 
The isophotes are disky within lh~^ kpc . The surface 
brightness profile of this source deviates from the R^ 
law. If we fit it with the i?^ law, we find a bump at 
1.5h~^ kpc and a dent at 2.8h~^ kpc. The inner bump 
is related to the inner arms. Apart from the outer part, 
this galaxy has a round appearance (e < 0.3). The rip- 
ples, loop structures and the deviation from the R^ law 
suggest that that IR 14378—3651 is a merging galaxy 
and is still relaxing towards the R^ light distribution. 



IR 10559+3845 From Fig. 10, IR 10559+3845 appears 
to have two nearly symmetric diffuse tails. NICMOS 
observations reveal two very close nuclei (Borne et al 
1998), indicating that this galaxy is a merging galaxy 
as well. The position angle of the major axis continu- 
ously changes by about 120 degrees from the center to 
the outer part; the position angle change in the inner 
part is more gradual than that in the outer part. The 
ellipticity of the inner isophotes is smaller (about 0.3) 
than the outer value (about 0.5). The center of the in- 
ner isophotes offsets to the south with respect to that 
of the outer ones. The mergering process appears not 
yet complete for this galaxy. 

IR 20176—4756 From both image and surface bright- 
ness profile (see Fig. 11), this galaxy appears to have 
been well relaxed. Remarkably, the surface brightness 
profile can be well fitted by a single exponential, rang- 
ing from 0.8/i^^kpc to 9.5ft.^^kpc. The central part 
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EPS files of the images are available at 
http://vega.bac.pku.edu.cn/~zhengz/preprint.html 
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Fig. 9. The left panel shows the HST Wide Field Camera image for IR 14378-3651 (the object labeled as GO). The 
image shown is 74.5" on a side. A satellite galaxy (Gl) and a foregromid star SI are visible. The middle panel shows 
the contours centered on IR 14378-3651; the contour levels are 1.9, 3.8, 7.6, 15.2, 20, 30.4, 40, 60.8, 70, 121.6, 243.2, 
486.4, 972.8, 1943.6 ADU, respectively. The thick horizontal bar indicates a scale of 5/i^^kpc. The north and east 
directions are indicated at the top right in the middle panel. The right panels show the variations of the ellipticity, 
position angle, surface brightness, and i?4/a (see §3.3) as a function oi R* . 

IR 10559+3845 



EPS files of the images are available at 
http://vega.bac.pku.edu.cn/~zhengz/preprint.html 
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Fig. 10. Same as in Fig. 9, for IR 10559+3845. The contour levels are 1.7, 3.4, 6.8, 13.6, 27.2, 54.4, 108.8, 217.6, 
435.2 ADU, respectively. 



(< lAh^^ kpc ) has boxy isophotes. On the other hand, 
there still exist merger signatures, such as a loop struc- 
ture (or inner arms) in the inner region, fan-like pro- 
trusion and thick tail-like material at the outer region 
(marked with an arrow in the image). The most signif- 
icant isophotal feature of IR 20176—4756 is that the 
position angle gradually decreases by about 180°. We 
speculate that the central boxy isophotes and isophotal 
twists may be caused by a central bar. 
IR 09039+0503 This galaxy has a remarkable mor- 
phology as shown in Fig. 12. A thick and luminous 
tail emerges from the north and bends to the south 
with a total length of about ~ 15/i~^ kpc . Two short 
tails locate to the south and northwest, respectively. 
A closer look reveals that this galaxy has two very 



close nuclei, which gives rise to the boxy isophotes 
within 1.2/i~^kpc. Although the central parts of the 
two galaxies have almost merged together, the galaxy 
is still far from complete relaxation and the surface 
brightness profile deviates significantly from the R^ 
law. There are many small and faint objects within 
30" of the source. It would be interesting to determine 
the redshifts of these small galaxies to see whether they 
are physically associated with the source. 
IR 23242—0357 It has a surface brightness profile sim- 
ilar to that of the second class, it has an inner com- 
ponent satisfying i?4 law and a more extended outer 
component. The primary difference is that this galaxy 
has a fainter nucleus. The inner region has two fan 
structures which cause the boxy isophotes at radius 
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EPS files of the images are available at 
http;//vega. bac.pku.edu.cn/~zhengz/preprint. html 
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Fig. 11. Same as in Fig. 9, for IR 20176-4756. The contour levels are 1.75, 3.5, 7, 9, 14, 20, 28, 35, 45, 56, 90, 112 
ADU, respectively. An excellent exponential fit to the surface brightness profile is also plotted. The arrow marks the 
protrusion and and thick tail-like material in the outer region. 



IR 09039+0503 



EPS files of tfie images are avaiiabie at 
http://vega.bac.pltu.edu.cn/~zhengz/preprint.html 
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Fig. 12. Same as in Fig. 9, for IR 09039+0503. The contour levels are 1.77, 3.54, 7.08, 14.16, 28.32, 56.64, 113.28, 
200, 226.56 ADU, respectively. An exponential fit to the surface brightness profile is also plotted. 



IR 23242-0357 



EPS files of the images are available at 
hLtp://vega. bac.pku.edu.cn/~zhengz/preprinL, html 
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Fig. 13. Same as in Fig. 9, for IR 23242-0357. The contour levels are 1.7, 3.4, 6.8, 13.6, 27.2, 54.4, 108.8, 217.6, 
435.2, 870.4 ADU, respectively. There are four other galaxies (G1-G4) within the field. The best R^ fit to the surface 
brightness profile is given by the straight line. 
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< 2.8/i~^kpc). The galaxy protrudes to the north- 
east and has a faint and sHm spiral tail around it in 
the south. There are also four fainter galaxies (Gl- 
G4) within the field, indicative of a group environment. 
Many fainter objects can be found within 15" of the 
source. 

To summarize, the five galaxies in class III have diverse 
appearances, including one (IR 14378—3651) with a close 
double nuclei. Their surface brightness profiles deviate to 
different degrees from the R* law, with IR 20176—4756 
fitted well by an exponential. IR 23242—0357 is similar to 
the objects in class II except that it has a fainter nucleus. 
We will return to the issue of the connection between the 
three classes of objects in the last section. 

5. The Optical Spectroscopic Properties 

Our studies so far have concentrated on the photometric 
properties of our sample galaxies. Clearly it is important 
to explore the spectroscopic properties of these galaxies in 
order to further understand these three classes of galaxies. 
In order to address this question and also to understand 
the possible group environment which we alluded to in 
Section 4, we have carried out spectroscopic observations 
from Dec. 20 to Dec. 23, 1998 using the 2.16m telescope 
at the Xinglong Station of Beijing Astronomical Observa- 
tory. We have observed IR 23140+0348, IR 09427+1929, 
IR 09039+0503 and IR 23242-0357, using a Zeiss uni- 
versal spectrograph mounted on the 2.16m telescope. A 
Tektronix 1024x1024 CCD was used giving a wavelength 
coverage of 4100A to 9100A with a grating of 200A/mm. 
The spectral resolutions are 9. 3 A (2 pixels). Wavelength 
calibration was carried out using a He-Ar lamp; the result- 
ing wavelength accuracy is better than lA. KPNO stan- 
dard stars were observed to perform flux calibrations. 

In the following section, we will present our observa- 
tional results together with data available in the literature, 
such as Kim et al (1998) and Lawrence et al (1999), for 
these three classes of objects. 

5.1. Class I 

The spectrum of IR 23140+0348 is shown in Fig. 14. 
As can be seen from this figure, there are strong emis- 
sion lines, as well as obvious 4000 A break (redshifted to 
~4900A), G-band and possible Call K and H absorption 
features. The Balmer absorptions are marginal in this low 
signal-to-noise ratio spectrum. From the emission line ra- 
tios and the classification criteria (Veilleux et al 1995; Os- 
terbrock et al 1989), IR 23140+0348 can be classified as 
a LINER. 

Compared with the "E+A"-like spectrum of ULIRGs 
FF J1614+3234 and TF J1020+6436 (Tran et al 1999; 
Breugel 1999; see also Dressier & Gunn 1983), the Balmer 
absorption line features of IR 23140+0348 is weaker. Nev- 
ertheless, the spectral features resemble each other and it 
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Fig. 14. The top panel shows the optical spectrum for IR 
23140+0348 . The bottom panel shows a magnified view 
for the spectral region between 4500A to 6000A. For com- 
parison of absorption lines, the spectrum of an A5 star is 
shown. Prominent lines are labelled. 



indicates that IR 23140+0348 has a population of young 
stars with age '^ 1 Gyr. Combining with the fact that 
its surface brightness profile is well described by the R^ 
law, we conclude that IR 23140+0348 is a young elliptical 
galaxy resulted from merging. This is further supported 
by the fact that IR 23140+0348 is a radio galaxy, and 
such galaxies are often associated with merger remnants 
(Barnes 1998). 

We have also observed three galaxies (GO, Gl and G2 
in Fig. 1) in the field of IR 09427+1929. We found that 
the redshifts for GO and G2 are 0.284 and 0.149, respec- 
tively; the low signal-to-noise ratio and the weak emission 
features make it difficult to securely identify Gl's redshift. 
The redshift of G2 is the same as the redshift given for IR 
09427+1929 by NEDQ The IRAS error ellipse (see Fig. 1) 
touches both GO and G2, with GO closer to the center of 
the error ellipse, so the identification is somewhat ambigu- 
ous. The spectrum, shown in Fig. 15, indicates that GO is 
a narrow line Seyfert 1 galaxy with extremely strong Fell 
emission. Since almost all extremely strong Fell emitters 
are ULIRGs, GO (at z = 0.284) is most probably the op- 
tical counterpart of the ULIRG. The NED identification 
seems to have confused the objects GO and G2. Incidently 



^ The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Admimistration. 
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Fig. 15. Spectrum for IR 09427+1929. Prominent lines Fig. 16. Spectrum for IR 23242-0357. Prominent lines 
are indicated. Notice the strong Fell lines. are labelled. 



G2 is a foreground galaxy at z = 0.149 with an Hll-region 
like spectrum. 

For the other two objects in class I, IR 04025—8303 
is also a Seyfert 1 galaxy with FWHM of Ha less than 
3000 km s~^ (Lawrence et al 1999) while the spectral clas- 
sification for IR 17432-5157 is LINER. 

5.2. Class II 

As for the spectra of galaxies in Class II, IR 00509+1225 
(I Zw 1) is a typical narrow line Seyfert 1 (NLSl) galaxy 
(Moran et al 1996); IR 10026+4347 is a NLSl-like QSO 
(Xia et al 1999); IR 00276-2859 is an IR QSO with 
strong Fell emission line and with FWHM of H/3 about 
2000 km s-i; IR 02055+0835 is also a QSO with strong 
Fell emission line and the FWHM of H/3 is less than 
3000 km s~^ from the QDOT redshit survey spectrum ob- 
tained by the William Hershel Telescope (Lawrence et al 
1999). It is remarkable that all four galaxies in Class II 
are classified as Seyfert 1 galaxies with strong Fell emis- 
sion lines and relatively narrow permitted emission lines. 
Furthermore, the optical luminosities for all four galaxies 
are in the range of QSOs. Except for IR 02055+0835 , all 
the other three objects were detected in the soft X-ray by 
ROSAT (the non-detection of IR 02055+0835 is perhaps 
because it is the most distant galaxy in the sample). All 
the three soft X-ray spectra are steep. Therefore, these 
four objects are NLSl-like QSOs (Boiler et al 1996), or 
they are en route to normal QSOs (Xia et al 1999). 



5.3. Class III 

The spectrum for IR 23242—0357 is a typical HII region 
spectrum (see Fig. 16). Table 2 gives the spectral classifi- 
cations for the other 4 galaxies in this class based on the 
data from the literature or the spectra from the QDOT 
redshift survey. IR 14378-3651 and IR 09039+0503 are 
classified as LINERs based on Kim et al (1998) and 
our spectrum (not shown), respectively. As Veilleux et 
al (1995) pointed out that the excitation mechanism for 
infrared luminous LINERs could be shock waves. There- 
fore it is possible that strong shocks are present in IR 
14378-3651 and IR 09039+0503, as in some starburst 
galaxies, such as NGC 6240 and NGC 3690 (Heckman et 
al 1990). To summarize, the statistical spectral properties 
of class HI objects are very different from those of Class 
I and Class II. These objects have either HII region-like 
or LINER-like spectrum, consistent with the presence of 
starbursts in these galaxies without obvious central AGNs. 

6. Summary and Discussions 

From a total of 58 HST WFPC2 Lband snapshot images of 
ULIRGs, we have selected 13 sources that are dominated 
by a single nucleus to study their surface photometry. Ac- 
cording to their surface brightness profiles, we find that 
these galaxies can be divided into three classes. Below we 
first briefly summarize their morphology, surface bright- 
ness profiles and optical properties, and then discuss the 
origin and evolution tracks of these galaxies. 

The surface brightness profiles of the four galaxies in 
class I are well fitted by the R^ law out to ~ lOft-^^ kpc . 
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with bumps at inner few kpc regions. The disky isophotes 
at the inner region <, l/i~^kpc are consistent with stel- 
lar disk formed by infalling gas during merging; molec- 
ular disks of comparable size have been detected in CO 
(e.g., Downes & Solomon 1998). The fine structures at the 
outer part of these galaxies hints that the relaxation pro- 
cess is incomplete and is slower than the inner part. The 
large position angle changes also hint that these galaxies 
are merger remnants and are still evolving into ellipticals. 
Each member of this class has round appearance. The 
round appearance suggests that they have experienced 
multiple mergers, since binary merging tends to produce 
more elongated objects (Weil & Hcrnquist 1996). These 
objects are at the last stage of merging. 

The four galaxies in Class II are composed of an inner 
component described by the R^ law out to ^ 5h^^ kpc 
plus some outer extension. The existing outer extension in- 
dicates that this class of ULIRGs is at intcrmidiate merg- 
ing stage and the merger remnant has not relaxed com- 
pletely. On the other hand, the outer extension may be 
from incompletely relaxed material or they can be stars 
formed in secondary infall from tidally ejected material 
(e.g., Barnes 1988). 

The five galaxies in class III show various degress of 
deviation from the R^ law, indeed IR 20176—4756 is very 
well fitted by an exponential law. Except IR 20176—4756, 
the other four ULIRGs are consistent with their being at 
an early merging stage, both from their morphologies and 
surface brightness profiles (see §4.3). 

The photometric properties and the optical spectro- 
scopic properties for our ULIRGs seem to be correlated. 
It is interesting that the percentage of Seyfcrt Is is very 
high (6/13) in our sample compared with complete ULIRG 
samples, such as Lawrence et al (1999), Kim ct al (1998), 
which has a percentage of Seyfert Is less than 15%. All 
class II galaxies are Seyfert 1 galaxies/QSOs while two 
out of four class I galaxies are Seyfert 1/QSOs as well. 
At least 5 of these 6 Seyfert Is /or QSOs are strong or 
extremely strong Fell emitters and most are NLSl like 
galaxies. This is very intriguing since our sample is solely 
selected using the criterion that there should be a single 
nucleus from HST snapshot images. In contrast, none of 
the Class III objects are Seyfert Is and they are classified 
as either HIIs or LINERs, which suggests that massive 
starbursts are present in these galaxies. 

Let us now investigate the evolutionary tracks between 
these three classes of ULIRGs. Morphologically speaking, 
class I objects are clearly the most relaxed, while class 
III objects in general show the most disturbed morpholo- 
gies and photometries, with class II objects somewhere in 
between. This suggests an evolution sequence from class 
III to class II, and then to class I, as relaxation proceeds 
to ultimately form a relaxed elliptical galaxy or SO galax- 
ies, depending on the initial mass ratios of the merging 
progenitors (Barnes 1998). 



In these three classes of galaxies, the interplay between 
star formation and central AGN activities is undoubtedly 
very complex, nevertheless, here we speculate qualitatively 
on the evolution of these galaxies based on such an in- 
terplay. First, as galaxies begin to interact and merge, 
massive nuclear starbursts are induced. In the mean time, 
some gas also flows into the center to ignite the central 
AGN activities, as indicated by numerical simulations. 
The spectroscopic properties of galaxies at this stage are 
a composite of the central AGN and circumnuclear star- 
bursts, with the latter being dominating. Our class III 
objects are probably at this stage. As time goes on, gas 
continuously fall into the center to fuel more intense cen- 
tral AGN activities, while the star formation has subsided 
due to the gradual depletion of the gas supply. Most of the 
dust is blown away, as a result, these galaxies are mostly 
seen as Seyfert 1 galaxies. As discussed by Xia et al (1999), 
the Seyfert 1 galaxy with strong Fell emission lines and 
with steep soft X-ray photon index could be at such a 
transition stage, namely these objects are emerging young 
QSOs. Our class II fits this description qualitatively. Fi- 
nally, as relaxation proceeds even further, gas is mostly 
consumed by star formation or blown out by superwind 
(Heckman et al 1990), and the central AGN weakens and 
the merger remnant emerges as a young elliptical galaxy 
- our class I objects. If the central gas supply is not yet 
exhausted, then some of these galaxies can also be Seyfert 
1 galaxies (as two of our class I objects are). This pic- 
ture was first pointed out by Sanders et al (1988a) ten 
years ago. Our observations seem to confirm this scenario. 
Along this sequence, the gas behavior is also expected to 
evolve. There is strong evidence that there are large (ten 
kpc scale), extended HI or emission-line nebula around 
some ULIRGs (Armus ct al 1990; Hibbard & Yun, 1997). 
The extended gas distribution is formed by the secondary 
infall of gas thrown out during the merging process. The 
infalling gas may form some outer extension of ULIRGs, 
as seen in our Class II objects. It will be very interesting 
to obtain HI information to study the evolution of gas in 
these galaxies (cf. Hibbard & Van Gorkom 1996). 

To conclude, from our analysis of the photometry of 
ULIRGs, we found that these galaxies are excellent lab- 
oratories and the surface brightness profile analysis is an 
important way to study the merging proccess. The obser- 
vational data provide a hint that some of these galaxies 
reside in groups of galaxies and multi-merger is not rare. 
Six out of these thirteen galaxies are Seyfert Is or QSOs, 
strongly suggesting that the AGN phenomenon is univer- 
sal in the merging process and is an integral part of el- 
liptical galaxy formation. A number of important works 
remain to be done. A multi-color campaign of these ob- 
jects will provide invaluable information about the star 
formation sites and the relative importance of starbursts 
and central nulcear activities. Furthermore, many galaxies 
(including a number of small and faint galaxies) are within 
the field of the primary ULIRGs. Whether these galaxies 
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are physically associated with the ULIRGs awaits future 
spectroscopic observations. Such observations will provide 
insight not only on the environments of ULIRGs but also 
on the nature of dwarf galaxies around these ULIRGs. 
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